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Abstract 

A jump-diffusion process is proposed to describe the dispiacements performed by singie myosin 
heads along actin filaments during the rising phases. The process consists of the superposition of 
a Wiener and a jump process, with jumps originated by sequences of Poisson-distributed energy- 
supplying pulses. In a previous paper 0, the amplitude of the jumps was described by a mixture of 
two Gaussian distributions. To embody the effects of ATP hydrolysis 1], we now refine such a model 
by assuming that the jumps' amplitude is described by a mixture of three Gaussian distributions. 
This model has been inspired by the experimental data of T. Yanagida and his co-workers concerning 
observations at single molecule processes level as described in |1U| and in the references therein. 

Keywords: Brownian motion, jumps, acto-myosin dynamics. 

A MS Classification: 92C17, 60J70. 

1 Introduction 

Muscle fibers are composed by a great number of even smaller fibers, called myofibrils, arranged parallel to 
the muscle fiber's major axis. The myofibrils consist of two different filaments: thick filaments, composed 
by myosin molecules, and thin filaments, composed by actin molecules. Myosin molecules' extremity 
stick out from the filament as small heads. Muscle contraction occurs when myosin molecules slide along 
actin filaments, fuelled by the chemical energy originating from ATP hydrolysis. The traditional model 
explaining the mechanism of myosin movement is the so-called "lever-arm swinging model" in which the 
neck region of the myosin head swings to generate displacement (see Cooke 0]). The swing motion is 
coupled tightly to the ATP hydrolysis cycle i.e. the myosin molecule moves in a single forward step 
during each ATPase reaction. The size of myosin displacement is about 6nm. 

Using new techniques for manipulating single actin filaments, Kitamura et al. (lUj have obtained 
highly precise and reliable measurements of the displacements performed by single myosin heads along 
actin filaments during their rising phases. They have recorded myosin displacements of lOnm ~ 30nm. 
These values indicate that, during each biochemical cycle of ATP hydrolysis, the myosin head may 
interact several times with an actin filament, undergoing multiple steps. The steps occur randomly in 
time, mainly in a forward direction, roughly not more than 10% in the backward direction. The step-size 
is approximately 5.5nm, that identifies with the interval between adjacent monomers in one strand of the 
actin filament. These observations clearly contradict the traditional lever-arm model, suggesting that a 
biased Brownian ratchet mechanism could be invoked: myosin head moves along the actin filament driven 
by the Brownian motion and the ATP hydrolysis biases the direction of the movement. 

In this paper we propose a stochastic model for the description of the displacements performed by 
the myosin head durig a rising phase, on account of the experimental results presented in Kitamura and 
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Yanagida [§]. Our model is not exactly a ratchet-based model but all the relevant basic ingredients are 
preserved: exploitation of Brownian motion, the assumption that energy released by ATP hydrolysis 
is responsible for the myosin displacements and the intrinsic asymmetry of the actomyosin system. To 
achieve directional motion, in addition to thermal noise and anisotropy, we consider an energy supply, 
resulting from the hydrolysis of ATP (see Buonocore and Ricciardi P]). We assume that this energy is 
stored by the myosin head and a part of it is released by small quanta of constant magnitudes. Usually 
not each energy release produces a step, but several releases are needed to generate a single step. We 
further assume that myosin head is embedded in a viscous fluid in thermal equilibrium and denote by (3 V 
the viscous friction drag coefficient that in the experimental conditions is estimated as 90 pN. 

Aiming to a phenomenological model to account for the available data on the displacements of the 
myosin head along an actin filament, we take as starting point the articles by Di Crescenzo et al. [S] 
and by Buonocore et al. in which a stochastic process consisting of the superposition of Wiener 
and jump processes is discussed. Myosin head's slide along the actin filament is viewed as a Brownian 
motion perturbed by jumps that occurr according to a Poisson process. Their amplitude is described by 
a mixture of three Gaussian random variables. 

The transition density of the stochastic process describing the stepwise motion and some of its mo- 
ments are obtained in Section [21 In Section we analyze the duration of the rising phase U and the 
position V of the myosin at the end of the rising phase. For both variables we obtain the density and 
some moments. In Section 0] the special case when the backward jumps are not allowed is considered. 
Our results appear to be in qualitative agreement with the available experimental observations. 



Let {X(t), t > 0} denote the real- valued stochastic process describing the motion of the myosin head 
along the actin filament during the rising phase. In the absence of jumps, the myosin head is assumed 
to move according to a Wiener process with zero drift and infinitesimal variance S 2 — 0.09 (nm) 2 /ms. 
Here use of relation 5 2 — 2kBT/f3 v has been made, with f3 v = 90 pN and where fcgT ~ 4pNnm is the 
thermal energy at absolute temperature T, ks = 0.01381 pN-nm/K denoting the Boltzmann constant and 
T = 293 K the environmental temperature. According to the previous assumptions, at the occurence of 
the i-th event of a Poisson process {N(i), t > 0} of intensity A, an energy release occurs and consequently 
the myosin head performs a jump of random amplitude Y{. The mean amplitude of such a jump is L, 
and — L with probabilities p, 1 — p — q and q, respectively. The following stochastic equation thus holds: 



where {B(t), t > 0} is the standard Brownian motion, and X(0) — 0. (We have arbitrarily set as the 
starting point of the motion). Note that the diffusive component of the motion has zero drift, according 
to the experimental evidence (see Kitamura et al. P0|)- We suppose that processes B(t), N(t) and r.v.'s 
Yi,Y2, • • ■ are independent, where YJ's are copies of a random variable Y characterized by distribution 
function Fy{y) and probability density The latter is assumed to be a mixture of three Gaussian 

densities fz 1 (y), fz 2 (y) and fz 3 (y) with means L, and — L respectively, and identical variance a 1 . 
Hence, for all y G R we have 
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with E(y) = L (p — q) and Var(F) = cr 2 . We set L = 5.5 nm (the distance between adjacent actin 
monomers) and a 1 = 0.1 (nm) 2 /ms. Densities /zi(y), fz 2 (y) an( A fz 3 (y) account for the following: 
fzt (y) describes the size of forward steps, fz 3 (y) that of backward steps, whereas fz 2 (y) accounts for an 
energy release that does not produce a step. 

For all x £ R and t > let us now introduce the following probability densities: 

f(x,t):=—P{X(t)<x\X(0) = 0}, (3) 
f n (x,t):=-^P{X(t)<x\X(0)=0,N(t)=n}, n=l,2,.... (4) 



One has: 



where 



f(x,t) = J2f n (x,t)P{N(t)=n}, (5) 



71 = 



P{N(t) = n} = ^4 . n = 0,l,..., (6) 



71 



is the probability distribution of a Poisson process N(t) having intensity A. The conditional density 
f n {x,t) can be expressed as follows: 

f w (x,t\x ) if n = 

fn{x,t) = { r00 (7) 

f w (x-y,t\x )dF^\y) tin = 1,2,..., 



where 

/w(M) = vi^ exp {~^}' xeR 

is the transition density of a Wiener process with zero drift and infinitesimal variance S 2 , and where 
F^"\y) denotes the n-fold convolution of Fyiy) with itself. From (J2J, by induction we obtain: 



d4' l) (V) = E ( J 0- - P - E ( J P 1 d ^ ( 8 ) 

fc=o ^ ' i=o 

where ipj,k,n{y) denotes a normal density with mean (2j — k)L and variance no 1 . Hence, from Q and 
®, for x G K, f > and n = 0, 1, ... we have: 



X 



Recalling Eq. JSJ, from Eqs. © and ©, for all a; g R and t > we finally obtain the probability 
density of X(t); 
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Note that, as shown in Di Crescenzo et al. jS], this is the solution of the following integro-differential 
equation: 

d f 5 2 d 2 f f°° 

with initial condition lim f(x,t) — S(x). Density l|l(J|l is multimodal, with peaks located at multiples of 

L (see the example plotted in Figure If p > q, at each point x = kL (k integer) where a peak exists, 
this peak is higher than the symmetric peak located at —kL, a manifestation of the prevalence of forward 
with respect to backward displacements. 




Figure 1: Plot of density JTOJ with a 2 = 0.1, p = 0.18, q = 0.01 and A = 1.10. 

Let us now obtain some moments of X(t). From we have: 

E[X(t)]:=m(t)=L(p-q)Xt, (11) 
Var[X(t)] :=v 2 {t) = {[L 2 (p + q)+a 2 ]\ + 5 2 }t, (12) 
E{[X(t) - m(t)] 3 } = L{p - q){L 2 + 3a 2 ) Xt. (13) 

Hence, from (I1H and l|12|) the coefficient of variation follows: 



cvwoi-^-^ y^A (14, 

m(t) L(p — qj\t 

Let us denote by fi(t) = X(t)/t the velocity of the myosin head displacement. Making use of Eqs. dllfl 
and JEJ we obtain the mean and the variance of nit): 

E[»(t)] = L(p-q)\, (15) 

v r , A} [L 2 (p + q)+a 2 }X + 6 2 
Var[^(*)] = . 

Finally, Eqs. (|T2^ and ifHfl) yield the skewness: 

E[X(t) - m(t)} 3 = X{p + g)L{p-g){L 2 + 3(T 2 ) 

v 3 (t) {{L 2 (p + q) + a 2 ]X + 5 2 } 3 / 2 y/t' 

The skewness goes to zero as t — * +oo, it is positive (negative) as p > q (p < q), whereas it vanishes 
when p — q. Under the assumption p > q, which is of interest in the context of myosin head's motion, 
we have the following monotonicity properties of the skewness: 

(i) as a function of A S (0, +oo), it is increasing for A < A = 2<5 2 /[L 2 (p + q) +a 2 }, decreasing for A > A, 
and goes to as A — > +oo; 



3 RISING PHASE 



5 



(ii) as a function of a 2 S (0, +oo), it is increasing if a 2 < a 2 = L 2 [2 (p + q) — 1] + 2<5 2 /A, decreasing for 
a 2 > cx 2 , and tends to when cr 2 — > +oo. 

Let us denote by M the r.v. describing the number of energy quanta releases required to produce one 
step, that corresponds to a jump of mean amplitude — L or L. On the ground of our assumptions, M is 
a geometric r.v. with parameter p + q. Denoting by D the r.v. describing the dwell time, i.e. the time 
between consecutive myosin jumps, we have 

D = T 1 + ---+T M , (17) 

where Ti is the duration of the effect of the z-th energy quantum release. Being N(t) a Poisson process 
with intensity A, we have that T\,T%, . . . are i.i.d. exponential r.v.'s with mean A . This is in agreement 
with Kitamura and Yanagida j^j, whose dwell-time histograms are well fitted by exponential curves. From 
(|17|l we obtain 

E(e^) = £ E(e^+-+^) P(M = *) = A( P + g) , s < X(p + q), 
so that the dwell time D is exponentially distributed with mean 

E(D, = W < 18 > 

From (|15|l and l|18(l we have 

p + q E{D) 

so that the mean velocity E[/x(t)] can be seen as the mean net displacement in one step divided by the 
mean time needed to have one step. 



3 Rising phase 

Some histograms of the number of steps per displacement have been shown in Kitamura and Yanagida 
[3] and Kitamura et al. |l(Jj . There are 66 observed displacements under a near-zero load and 77 
displacements for a load between and 0.5 pN. The number of net steps in each displacement ranges 
from 1 to 5 under a near-zero load and from 1 to 4 for a load between and 0.5 pN. Moreover, the 
distribution of the total number of net steps observed in a rising phase minus 1 is well-fitted by a Poisson 
random variable. Hence, in order to include such features in our model we introduce a random variable 
U that describes the duration of the myosin head rising phase. Let {N + (t), t > 0} and {N-(t), t > 0} 
be two independent Poisson processes characterized by intensities Xp and Xq, that describe the numbers 
of forward and backward steps performed by the myosin head during a rising phase, respectively. We 
assume that U is the first-passage time of the number of net steps, given by R(t) := N+(t) — N-(i), 
through a Poisson-distributed random threshold S, with 

P(S=k) = G P ? ] , fc = l,2 (19) 

where p > (note that this assumption refines a previous model considered in Buonocore et al. 2). The 
mean number of net steps during the rising phase is thus given by 

E(S)=p + l. (20) 

We note that R(t) is a randomized random walk, whose transition probability for positive A, p and q is 
given by (see, for instance, Conolly 

Pk(t) = P{R(t) = k | R(0) = 0} = ( - ) e- x ^+^ h^Xty/pql t>0, 
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where Ik{x) — ^2^=o {x/2) 2n+k /[n!(n + fc)!] is the modified Bessel function. For all t > 0, the p.d.f of U 
can be expressed as follows: 

+00 

fu(t) = Y i fu\s(t\k)P(S = k), (21) 

k—l 

where P(S — k) is given in (|T!)|) and fu\s{u | k) is the first-passage-time density of R(t) through k. The 
latter is given by 



fu\s(t\k) = * Pk (t) = * (?) ' e- x ^ I k (2\t^j). t>0, 
From (Unj, (EU and we have 



fu(t) 



-\{p+q)t +2^ £ 



k/2 



'^Jj—TyXl) p^^Vm), t>o. 



Making use of Eq. lt^|) and of Eqs. 4.16.1 and 4.16.2 of Erdelyi et al. [7], for p > q we obtain: 



E(U) = 



Var(ZJ) = 



CV(U) = 



X(p-qY 
(2 P +l)p + q 



(22) 



(23) 



(24) 



A 2 (p — q) 3 



1 (2p+l);p + g 



P+ 1 V P~ 9 




Figure 2: Plots of density IL'MI in the following cases: (i) p = 1.5, [X(p + q)] 1 = 4.8 ms and p/q — 18, where 
E(U) = 13.37ms and Var(C7) = 122.80ms 2 , and (ii) p = 0.6, [\(p + q)]' 1 = 12.2ms and p/q = 6, where 
E(U) = 27.83 ms and Var((7) = 858.98 ms 2 . 

Note that E(C/) = E(S)/Fj[R(t)/t], so that the mean duration of the rising phase equals the mean number 
of net steps times the net steps rate. Two plots of density (|23|) are shown in Figure 2. Parameters p and 
q are chosen according to the experimental results concerning myosin number of steps along the actin 
filament for different load values that are given in Kitamura and Yanagida [2], where the ratio between 
the number of forward steps and backward steps is 18 under a near-zero load, and is 6 for a load between 
and 0.5 pN. Moreover, recalling from (|18fl that the mean dwell-time is [X(p + q)]~ , A is obtained from 



[Xip + q)}-^ 



12.2 ms 



under a near-zero load 

for a load between and 0.5 pN 
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where the rates of dwell times, appearing on the right-hand sides, are experimentally evaluated in the 
presence of 1 pM ATP and at 20 C (see Kitamura et al. JO] an d Kitamura and Yanagida [S]). Moreover, 
due to (|20|) . the value of p is chosen according to the following experimental mean numbers of steps 
observed during rising phases (see Kitamura and Yanagida [§]): 



2.5 under a near-zero load 

1.6 with a load between and 0.5 pN. 



Hence, in case (i) of Figure 2, where fu{t) exhibits a positive mode, it is p/q — 18, X(p + q) = (4.8) 1 
and p = 1.5; in case (ii), where fu(t) is strictly decreasing, we have p/q — 6, A(p + q) = (12.2) -1 and 
/; = 0.6. 

Let us now denote by V the random variable denoting the position attained by the myosin head at 
the end of the rising phase. Its density is given by: 



fv(x) 



f(x,t)fu(t)dt, 



x e 



(25) 



with f(x,t) expressed in (|10[1 and fu(t) given in (12,'SH . In order to provide a qualitative insight of its 
behaviour, in Figures 3 and 4 we show two instances of probability density l|25|l . For p > q, making 
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Figure 3: Plot of density ^ with p = 1.5, a 2 = 0.1, p = 0.72, q = 0.04 and A = 0.27. 
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Figure 4: Plot of density ^ with p = 0.6, a 2 = 0.1, p = 0.18, q = 0.03, A = 0.39. 

use of Eq. H25|l . we have: 

E(V) = (p+l)L, 



Var(V") = 



XL 2 [{2p + l)p + q] + (p + l){A[L 2 (p + q) + a 2 } + S 2 } 
X(p-q) 



(26) 
(27) 
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Table 1: See text for explanation 



p 


q 


A 


E(V) - 2 v / Va,r(V) 


E(V) + 2 n /Var(V) 


0.09 


0.005 


2.19 


-15.82 


43.32 


0.18 


0.01 


1.10 


-15.72 


43.22 


0.36 


0.02 


0.55 


-15.67 


43.17 


0.72 


0.04 


0.27 


-15.65 


43.15 



CY(V) = y/AL^ [(2p + \)p + q] + (p + l){\[L*(p + q) + g] + g} 

Cp + l)VA(p-g) 

Comparing Eq. f2*6^l with l|2*4*jl and lfT5| we have E(V) = E(t/) ■ E[/t(i)], so that the mean position at the 
end of the rising phase equals the mean duration of the rising phase times the constant mean velocity of 
the myosin head. In Tables ^and [21 the values of E(V) ±2-y/Var(V) are indicated. These are obtained via 
Eqs. it2*6")l and lt2*Tjl for some choices olp, q and A under the conditions p/q = 18, p = 1.5 and E(V) = 13.75 
(Table [TJ and p/q = 6, p = 0.6 and E(V) = 8.8 (Table HJ. 



Table 2: See text for explanation 



p 


q 


A 


E(V) - 2y/Va,T{V) 


E(V) + 2y/Var(V) 


0.03 


0.005 


2.34 


-16.70 


34.30 


0.06 


0.01 


1.17 


-16.45 


34.05 


0.12 


0.02 


0.59 


-16.32 


33.92 


0.24 


0.04 


0.29 


-16.26 


33.86 



4 Rising phase in the absence of backward steps 

The experimental results of Kitamura and Yanagida jH] and Kitamura et ol. |10| show that backward 
steps during a rising phase are rare with respect to forward steps. Hence, it is challenging to analyse 
the results given in the previous section when backward steps are not allowed, i.e. by assuming q = 0. 
In this section we thus assume that R(t) identifies with iV+(t), which is a Poisson process with intensity 
A p. Hence, the random duration U of the rising phase in the absence of backward steps is the first- 
passage time of N+(t) through the random threshold S. Recalling Eq. I)18[l. in this case the dwell time is 
exponentially distributed with mean value (Ap) _1 , so that fu\s(t I k) is an Erlang density with parameters 
k and (Ap) _1 : 

f uls (t\k) = Xpe-^^^-, t>0. (28) 

From Eqs. (Tffi , (|2~T|) and (J2BJ, we obtain the probability density of the rising phase duration in the 
absence of backward steps: 

fu (t) = e~" Xp e- Xpt I (2 Xp t) , * > 0. (29) 

From l|29|l we have that /c/(0) = e~ p Xp > 0, and that fu{t) is a Polya frequency of order 2 density (PF2), 
i.e. a logconcave density (see Marshall and Olkin JI]). This means that [U — t\ U > t] >i r [U — r | U > r] 
whenever < t < t (see Theorem 1.C.22 of Shaked and Shantikumar ^2]). In other words, as time goes 
on, the residual time of rising phase decreases in the likelihood ratio order sense. Moreover, from l|29|) we 
obtain 

= c-> (\p) 2 e-W [(p-l)I (2y/p~\Ji) - P I 2 {2^p~XVtj 



t > 0. 
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From this expression it is possible to show that density fu(t) is decreasing for all t > if p < 1, whereas 
density (|29|) is unimodal with a positive mode if p > 1. 

Acknowledgments 

Work performed within a joint cooperation agreement between Japan Science and Technology Corporation 
(JST) and Universita di Napoli Federico II, under partial support by INdAM (G.N.C.S.). 

References 

[1] A. Buonocore and L.M. Ricciardi. Exploiting thermal noise for an efficient actomyosin sliding mech- 
anism. Mathematical Biosciences, to appear, 2003. 

[2] A. Buonocore, A. Di Crescenzo and B. Martinucci. A stochastic model for the rising phase of myosin 
head displacements along actin filaments. Preprint, 2002, Dipartimento di Matematica e Appli- 
cazioni, University of Naples Federico II. 

[3] B. Conolly. "Lecture Notes on Queueing Systems." Ellis Horwood, Chichester, 1975. 

[4] R. Cooke. Actomyosin interaction in striated muscle. Physiological Reviews 77: 671-697, 1997. 

[5] D. Cyranoski. Swimming against the tide. Nature 408: 764-766, 1999. 

[6] A. Di Crescenzo, B. Martinucci and L.M. Ricciardi. On modeling the rising phase of myosin head 
displacements in single molecules processes. Cybernetics and Systems 2002, 16th European Meeting 
on Cybernetics and Systems Research, Austrian Society for Cybernetics Studies, Vienna: 295-300, 
2002. 

[7] A. Erdelyi, W. Magnus, F. Oberhettinger and F.G. Tricomi. "Tables od Integral Transforms." Vol. 
I. McGraw-Hill, New York, 1954. 

[8] Y. Ishii and T. Yanagida. Single molecule detection in life science. Single Molecules 1: S.5-16, 2000. 

[9] K. Kitamura and T. Yanagida. Stochastic properties of actomyosin motor. BioSystems (submitted). 

[10] K. Kitamura, M. Tokunaga, A.H. Iwane and T. Yanagida. A single myosin head moves along an 
actin filament with regular steps of 5.3 nanometres. Nature 397: 129-134, 1999. 

[11] A.W. Marshall and I. Olkin. "Inequalities: theory of majorization and its applications." Mathematics 
in Science and Engineering, Vol. 143. Academic Press, New York, 1979. 

[12] M. Shaked and J.G. Shantikumar. "Stochastic Orders and Their Applications." Academic Press, San 
Diego, 1994. 

A. Buonocore, B. Martinucci and L.M. Ricciardi: DIPARTIMENTO DI MATEMATICA E APPLICAZIONI, 
Universita di Napoli Federico II, Via Cintia, Napoli 1-80126, Italy 

A. Di Crescenzo: Dipartimento di Matematica e Informatica, Universita di Salerno, Via S. 
Allende, 1-84081 Baronissi (SA), Italy 



